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We examine a method to determine the neutron-skin thickness of nuclei using data on the charge-exchange anti-analog giant dipole 
resonance (AGDR). Calculations performed using the relativistic proton-neutron quasiparticle random-phase approximation (pn- 
RQRPA) reproduce the isotopic trend of the excitation energies of the AGDR, as well as that of the spin-flip giant dipole resonances 
(IVSGDR), in comparison to available data for the even-even isotopes '^'^Sn. It is shown that the excitation energies of the 
AGDR, obtained using a set of density-dependent effective interactions which span a range of the symmetry energy at saturation 
density, supplemented with the experimental values, provide a stringent constraint on value of the neutron-skin thickness. For 
'^^^Sn, in particular, we determine the value A/?p„=0.21+0.05 fm. The result of the present study shows that a measurement of the 
excitation energy of the AGDR in {p,n) reactions using rare-isotope beams in inverse kinematics, provides a valuable method for 
the determination of neutron-skin thickness in exotic nuclei. 
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1. Introduction 

■ An interesting phenomenon in nuclear structure is the forma- 
tion of a skin of neutrons on the surface of a nucleus, and its 
evolution with mass number in an isotopic chain [IJ. A precise 
^measurement of the thickness of neutron skin is important not 
only because this quantity represents a basic nuclear property, 
,but also because its value constrains the symmetry energy term 
of the nuclear equation of state A detailed 

knowledge of the symmetry energy is essential for describing 
the structure of neutron-rich nuclei, and for modeling proper- 
ties of neutron-rich matter in applications relevant for nuclear 
astrophysics. 

', The diff'erence between the neutron and proton rms radii 
is rather small (few percent) and a precise measurement of 
the neutron-skin thickness presents a considerable challenge. 
■Several methods have been used to determine this quantity 
but almost all of these are applica- 
ble only to stable nuclei and the results are model dependent 
ifioll . Methods based on coherent nuclear motion include ex- 
citations of the isovector giant dipole resonance (IVGDR) |@], 
the isovector spin giant dipole resonance (IVSGDR) [W], the 
Gamow-Teller resonance (GTR) measured relative to the Iso- 
baric Analog State (IAS) lll4l|. and high-resolution study of the 
electric dipole polarizability |3tl- 

The Pb Radius Experiment (PREX) at JLAB [5] has initi- 
ated a new line of research based on the parity-violating elastic 
electron scattering to measure the neutron density radius R„, 



which in turn allows to determine the neutron-skin thickness 
from Mi 



pn 
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R„ - Rp, where Rp is the radius of the proton 
density distribution. Although parity-violating elastic electron 
scattering provides a model independent measurement of ARp„, 
its current precision is far from satisfactory and the method can- 
not be applied to unstable isotopes. 

Radioactive ion beams (RIBs) have recently been employed 
to determine the neutron-skin thickness in unstable nuclei, 
specifically in measurements of reaction cross sections and 
pygmy dipole resonances For an accurate determi- 

nation of this quantity using RIBs, it is imperative to find a fea- 
sible method that employs reactions with low-intensity RIBs in 
inverse kinematics. We have recently introduced a new method 
1 1/71 uM based on the excitation of the anti-analog giant dipole 
resonance (AGDR) observed in {p,n) reaction lll9ll . As pointed 
out by Krmpotic, the excitation of the AGDR depends sensi- 
tively on the neutron-skin thickness f?^ and, therefore, ARpn 
could be deduced from the measurement of AGDR excitation 
energy. 

The main objective of this work is to test the method that de- 
termines the neutron-skin thickness in nuclei from AGDR data. 
By calculating excitation energies E(AGDR) and ARp„ in a fully 
self-consistent theoretical approach, and comparing to available 
data, the feasibility of the method will be tested with the aim to 
provide a basis for future studies with RIBs. 

The paper is organized as follows. In Sec. 2 we present a 
short outline of the AGDR mode and discuss its properties with 
the aim to constrain the neutron-skin thickness. The experi- 
mental results on AGDR and IVSGDR obtained from ip,n) and 
(^He,f) reactions for Sn isotopes are reviewed in Sec. 3. The 

February 26, 2013 



theoretical framework and model calculations of excitation en- 
ergies of the AGDR and IVSGDR, and neutron-skin thickness 
for the even Sn-isotopes, are described in Sec. 4. The centroid 
of the AGDR for '^"^Sn is deduced from an earlier {p,n) mea- 
surement in Sec. 5, and used to determine the corresponding 
neutron-skin thickness by comparing with model calculations. 
Section 6 includes the summary and a brief outlook for future 
studies. 

2. Charge-exchange AGDR and neutron-skin thickness 

The AGDR corresponds to the AJ" = 1 , AL = 1 reso- 
nant excitation, and represents the anti-analog giant dipole res- 
onance, i.e., the Tq-I component of the charge-exchange GDR 
(ro is the ground-state isospin of the target nucleus). Figure[T]il- 
lustrates the ground state and the giant dipole resonance (GDR) 
state of a target nucleus (T, = To), and the corresponding reso- 
nant states in the daughter nucleus reached by the ip,n) charge- 
exchange reaction: the IAS (isospin = To), and the anti-analog 
(isospin = To - 1) states: GTR, IVSGDR, and AGDR 1 19]. 

The transition strength of dipole excitations is fragmented 
into the Tq-I, To and To -i- 1 components because of the isovec- 
tor nature of the (p,n) reaction. The pertinent Clebsch-Gordan 
coefficients [2T] show that the To - 1 component (AGDR) is fa- 
vored with respect to the To and To + I components, by factors 
of To and 2T^, respectively. Accordingly, Fig. [1] shows only 
components that are strongly favored by isospin selection rules 
for large To 1.191 . 

For the charge-exchange dipole operator 



nYioifi)TT, 



(1) 



where t^I^ denotes isospin raising and lowering operators, one 
obtains the non-energy-weighted sum rule (NEWSR) |22], 



S- -S^ ^—(N{rl)-Z{rl)). 



(2) 



S " and S ^ denote the sums of transition strengths in the and 
channels calculated using Eq. ([T]i, respectively. The AGDR 
mode is mediated by the O operator We note that the sum 
rule Eq. (|2]l is proportional to the one of the IVSGDR (up to 
the factor three that corresponds to different spin components 
in the latter case}, previously used to determine the neutron- 
skin thickness lllOl 12311 . In Ref. ll22ll . Auerbach et al. derived 
an energy-weighted sum rule (EWSR) for the dipole strength 
excited in charge-exchange reactions. However, this EWSR is 



not sensitive to the neutron-skin thickness 02211 . 

Although the NEWSR provides information on the neutron- 
skin thickness, in practice it is not straightforward to deter- 
mine ARp„ from the sum rule. The problem is that usually the 
non-energy weighted strengths are not simultaneously available 
for both and channels, i.e., experiments are mainly per- 
formed in the channel and provide S " only. To deduce ARpn 
from the NEWSR additional approximations or theoretical in- 
put are needed, e.g. the estimate for 5^ and the normalization 
constant ifioll . One can, however, intuitively understand the re- 
lation between the neutron-skin thickness and the energy of the 



AGDR if one considers that for nuclei with - Z » 1 one 
can neglect S ^ because of Pauli blocking. In such a case, the 
EWSR, which is a constant, can be expressed as the product 
of the NEWSR and the AGDR energy and thus one can easily 
understand the inverse proportionality between the AGDR en- 
ergy and neutron-skin thickness, i.e. Eacdr decreases if ARpn 
increases and vice versa. 

Instead of the NEWSR constraint on ARp„, in this work we 
aim to establish an alternative approach motivated by the study 
of Ref. 120], where a simple schematic model indicated strong 
sensitivity of the AGDR excitation energy on ARp„. In Sec. 4 
we will explore this relation using a fully microscopic theoreti- 
cal approach and available data on the AGDR energies. 
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Figure 1: The ground state and the GDR of the target nucleus (7, = Tq). Also 
shown are the IAS (isospin=ro) and anti-analog states (isospin=ro - 1) GTR, 
IVSGDR and AGDR in the daughter nucleus (T- = Tq - I), excited in a {p,n) 
reaction. Only these components are shown because they are strongly favored 
by isospin selection rules for large Tq. See also Ref. 11911 . 



3. Isovector giant resonances excited by (p,n) reactions 

The first identification of a giant dipole transition excited in 
charge-exchange (p,n) reactions was reported by Bainum et al. 
[24] for the case ''°Zr(/?,n)''°Nb at 120 MeV In addition to the 
pronounced excitation of the GTR, a broad peak was observed 
at an excitation energy of 9 MeV above the GTR, with an an- 
gular distribution characteristic of a AL - 1 transfer This ex- 
citation energy is about 4 MeV below the location of the T = 5 
analog of the known GDR in ^"Zr, and thus it was suggested 
that this state is the T = 4 anti-analog of the GDR. 

Dipole resonances have also been studied systematically in 
(p,n) reactions at E,, = 45 MeV by Sterrenburg et al. lll9ll . 
for 17 different targets from ^^Zr to ^"^'Pb. Nishihara et al. 
lEsIl measured also the dipole strength distributions at E,, =41 
MeV. It was shown |2^ 27] that the observed AL- 1 resonance 
in general corresponds to a superposition of all possible spin- 
flip dipole (IVSGDR) and non-spin-flip dipole (AGDR) modes. 
According to Osterfeld [21], the non-spin-flip to spin-flip ratio 
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is favored at low bombarding energy (below 50 MeV), and also 
at very high bombarding energy (above 600 MeV). Properties 
of the IVSGDR were further investigated by Gaarde et al. [128 1 
using ip,n) reactions on targets with mass of 40 < A < 208, 
and by Pham et al. 1I29I1 in (^He,f) reactions. In every experi- 
mental spectrum a peak was observed at an energy several MeV 
above the GTR, with an angular distribution characteristic of a 
AL - 1 transfer. The observed excitation energies of the AGDR 
lfl9ll and IVSGDR Qlli], relative to the IAS, are shown in 
Fig. 2 as functions of the mass number for the even-even Sn 
isotopes. For both modes one observes a systematic decrease 
of the excitation energy along the Sn isotopic chain. 



Gamow-Teller gian t resonance was studied by the (p,n) reac- 
tion with RIBs 13311 using a similar neutron spectrometer ll34ll . 
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Figure 2: Excitation energies of the AGDR and IVSGDR relative to the IAS 
for the even-even Sn isotopes as functions of the mass number. The circles 
represent the experimental results from Sten'enburg et al. 1 19], the squares 
are from Pham et al. |29], the triangles from Krasznahorkay et al. | 23,], and 
the stars are the pn-RQRPA values for the AGDR (higher) and the IVSGDR 
(lower), calculated with the DD-ME2 effective interaction t30i1 . 

The {p,n) reaction has a high cross-section, and in inverse 
kinematics the energy of the neutrons is only a few MeV, which 
can be measured with highly efficient detectors. In our recent 
experiments (vf\, a 600 MeV/nucleon '^"^Sn relativistic heavy- 
ion beam was directed onto a hydrogen target. The ejected neu- 
trons were detected by a low-energy neutron-array (LENA) ToF 
spectrometer f3l. 32 1, developed in Debrecen. The spectrom- 
eter was placed at 1 m from the target, covering a laboratory 
scattering-angle region of 65° < &lab ^ 75°. In this way, 
the excitation energy of the AGDR and IAS can be determined, 
and new data enable studies of neutron-skin thickness in nuclei. 
One expects that future progress with RIBs and novel experi- 
mental techniques should provide data on the AGDR in exotic 
nuclei with even more pronounced neutron skin. The exper- 
imental feasibility of the suggested method is also supported 
by a recent publication in which the strength distribution of the 



4. Theoretical analysis 

To describe the evolution of excitation energies of the AGDR 
relative to the IAS, and their relation to ARp„, we perform a mi- 
croscopic theoretical analysis based on relativistic nuclear en- 
ergy density functionals. The theoretical framework is realized 
in terms of the fully self-consistent relativistic proton- neutron 
quasiparticle random-phase approximation (pn-RQRPA) based 
on the relativistic Hartree-Bogoliubov model (RHB) |35]. The 
pn-RQRPA is formulated in the canonical single-nucleon basis 
of the RHB model in Ref. |36], and extended to the descrip- 
tion of charge-exchange excitations (pn-RQRPA) in Ref. |37]. 
The RHB + pn-RQRPA model is fully self-consistent: in 
the particle-hole channel effective Lagrangians with density- 
dependent meson-nucleon couplings are employed, and pairing 
correlations are described by the pairing part of the finite- range 
Gogny interaction [38]. 

For the purpose of the present study we employ a family 
of density-dependent meson-exchange (DD-ME) effective in- 
teractions for which the constraint on the symmetry energy at 
saturation density was systematically varied, and the remain- 
ing model parameters were adjusted to reproduce empirical 
nuclear-matter properties (binding energy, saturation density, 
compression modulus), and the binding energies and charge 
radii of a standard set of spherical nuclei ll39ll . These effective 
interactions were used to provide a microscopic estimate of the 
nuclear-matter incompressibility and symmetry energy in rela- 
tivistic mean-field models i39(l . and in Ref. 1 15] to study a pos- 
sible correlation between the observed pygmy dipole strength 
in '^'^•'^^Sn and the corresponding values for the neutron-skin 
thickness. In addition to a set of effective interactions with sys- 
tematically varied values of the symmetry energy at saturation 
density, the relativistic functional DD-ME2 130(1 is also used 
here to calculate the excitation energies of the AGDR with re- 
spect to the IAS, as a function of the neutron skin. Pertinent 
to the present analysis is the fact that the relativistic RPA with 
the DD-ME2 effective interaction predicts for the dipole polar- 
izability ll 
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(directly proportional to the inverse energy-weighted moment 
m_i) of 208pb the value 0-0=20.8 fm^, in very good agreement 
with the recently measured value: = (20.1 + 0.6) fm-' [3j]. 

In addition to the experimental excitation energies. Fig. 2 
also includes the theoretical results obtained with the RHB -H 
pn-RQRPA model using the DD-ME2 effective interaction. The 
difference in the excitation energy of the AGDR and the IAS, 
as well as between the IVSGDR and the IAS, for the even-even 
Sn isotopes are shown as functions of the mass number For 
the excitation energies of the AGDR and the IVSGDR we take 
the centroids of the theoretical strength distributions: mi /mo, 
whereas a single peak is calculated for the IAS. Within the 
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experimental uncertainty, we find a reasonable agreement be- 
tween the data and the theoretical values for the AGDR. The 
largest deviations, x 0.4 MeV, correspond to '^^•'^'*Sn. The 
agreement is less satisfactory for the IVSGDR, with the dis- 
crepancy being especially large for ""*Sn and "*Sn. Part of 
the difference is probably caused by the overlapping of the two 
resonances. Experimentally it is possible to enhance the exci- 
tation of the AGDR with respect to the IVSGDR by choosing 
a low bombarding energy (< 45 MeV), but it is expected that 
the suppressed IVSGDR can still cause some lowering of the 
energy centroid of the AGDR observed in the (p,n) reaction 
since the former has a lower excitation energy. On the other 
hand, the excitation of the IVSGDR can be enhanced at higher 
bombarding energy (around 150 - 200 MeV) but a small frac- 
tion of the AGDR still remains, raising the centroid energy as 
the energy of the AGDR is higher than that of the IVSGDR. 
Since the IVSGDR has three components with AJ" - , 1~ 
and 2 , its strength generally spreads over a larger energy in- 
terval compared to the AGDR. The model calculation cannot 
reproduce these structures as precisely as the centroid energy 
of the AGDR, and this is the reason why in this work we make 
use of the AGDR to determine the neutron-skin thickness. 

The calculated values of the neutron-skin thickness for the Sn 
isotopes as a function of mass number are compared to avail- 
able data in Fig. 3. The RHB neutron-skin thicknesses ob- 
tained using the DD-ME2 effective interaction are in very good 
agreement with the experimental values obtained using vari- 
ous methods |2, JJJ. The self-consistent RHB calculation 
of ARp„, and the corresponding pn-RQRPA excitation energies 
of the AGDR, establish a connection between these quantities 
and suggest a feasible method for determining the neutron-skin 
thickness from AGDR data. 

5. Determination of the neutron-skin thickness of '^''Sn 

In this section the measured AGDR excitation energy for 
'^"^Sn, together with the consistent RHB plus pn-RQRPA model 
calculation of ARp„ and the AGDR energy, is used to constrain 
the value of the neutron-skin thickness. We consider the avail- 
able data for the AGDR for '^'^Sn from Sterrenburg et al. ^ 
(E(AGDR) - E(IAS) = 10.60 ± 0.20 MeV), but slightly in- 
creased to E(AGDR) - E(IAS ) = 10.93 + 0.20 MeV in order to 
approximately compensate the effect of the energy shift caused 
by the mixing with the IVSGDR. Below we explain how this 
energy shift is determined. 

Austin et al. |27] developed a phenomenological model to 
describe the variation with bombarding energy of the peak po- 
sitions of the AGDR and IVSGDR observed in (p,n) reactions. 
They assumed that the position C of the centroid of the AL = 1 
excitations (including both the AGDR and IVSGDR) at a bom- 
barding energy Ep is given by the weighted average of the en- 
ergies: 



C 



Q-qEq +O-1E1 
(To + CTi 



Eo- 



1 + CTi /(To 



A, 



(4) 




where Eo(Ei) is the energy of the AGDR (IVSGDR), A = £0 - 
£1 and (To((Ti) is the cross section for S = (S = 1) transfer 
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Figure 3: Calculated values of the neutron-skin thickness for the even-even Sn 
isotopes as a function of the mass number (filled circles connected by the solid 
line), compared to experimental results obtained with the antiproton absorption 
method 1 11] (triangles), from (p,p) scattering data | 2] (crosses), and with the 
IVSGDR method 110.1 normalized to the {p,p) result for '-''Sn {2^ (squares). 



They estimated the o-i/o-q ratio by : o-^/o-q ~ (Ep(MeV)/55)^ 
1 27|] and obtained the energy of the AGDR in ^^*Sn to be 14.4 
+ 2.2 MeV, which is completely different from any theoretical 
prediction 11271] . In reality, the centroid of the dipole strength 
distribution is usually determined by fitting the distribution by 
a Gaussian or a Lorentzian curve. This makes a significant dif- 
ference in case of '-'*Sn, where the AGDR and the IVSGDR 
display very diff^erent widths: 3.6 MeV HI and 9 MeV iS . 
respectively. 

To determine the energy shift of the AGDR peak at Ep = 45 
MeV from the empirical peak energy, we simulate the mixing 
of the AGDR and IVSGDR by using their real widths of 3.6 
MeV and 9 MeV, the ratio of their intensities as approximated 
by Austin et al. and their energy difference A=2.3 MeV 

obtained from Fig. 2. The composite spectrum is then fitted by 
a Gaussian curve in a reasonably wide energy range (+5 MeV) 
around the position of the peak, and this yields an energy shift 
of 0.33 MeV for the AGDR. 

The sensitivity of the centroid energy of the AGDR to the 
neutron-skin thickness of '^^Sn is explored by performing RHB 
-I- pn-RQRPA calculations using a set of the effective interac- 
tions with different values of the symmetry energy at satura- 
tion: 04 = 30, 32, 34, 36 and 38 MeV (and correspondingly 
different slopes of the symmetry energy [6j) and, in addition, 
the DD-ME2 effective interaction (04 = 32.3 MeV). In Fig. 4, 
the resulting energy differences E(AGDR) - E(IAS ) are plotted 
as a function of the corresponding neutron-skin thickness ARp„ 
predicted by these effective interactions. 

The two parallel solid lines in Fig. 4 delineate the region of 
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Figure 4: The difference in the excitation energy of the AGDR and the IAS for 
the target nucleus '-"^Sn, calculated with the pn-RQRPA using five relativistic 
effective interactions characterized by the symmetry energy at saturation 04 = 
30, 32, 34, 36 and 38 MeV (squares), and the interaction DD-ME2 (04 = 32.3 
MeV) (star). The theoretical values E(AGDR)-E(IAS ) are plotted as a function 
of the corresponding ground-state neutron-skin thickness ARp„, and compared 
to the experimental value. 



theoretical uncertainty for the used set of effective interactions. 
When adjusting the parameters of these interactions 113 Oi 13911 . 
an uncertainty of 10% was assumed for the difference between 
neutron and proton radii for the nuclei "^Sn, '^'*Sn, and ^*'**Pb. 
This set of interactions was also used to calculate the electric 
dipole polarizability and neutron-skin thickness of ^'"^Pb, '^^Sn 
and '^'^Ca, in comparison with the predictions of more than 40 
non-relativistic and relativistic mean-field effective interactions 



a. 

By comparing the experimental result for E{AGDR)-E(IAS ) 
to the theoretical energy differences (see Fig. 4), we deduce 
the value of the neutron-skin thickness 
0.205 ± 0.050 fm (including theoretical uncertainties). In Ta- 
ble I the value for AR„p determined in the present analysis is 
compared to previous results obtained with a variety of exper- 
imental methods. The very good agreement with previously 
determined values reinforces the expected reliability of the pro- 
posed method. 



6. Conclusion and outlook 

A method to determine the size of the neutron-skin thickness 
in nuclei using data on the anti-analog giant dipole resonance 
has been discussed. Charge-exchange (p,n) reactions provide 
an excellent probe for the neutron-skin thickness, as already 
demonstrated by measurement of the IVSGDR and GTR, and 
the AGDR provides a complementary approach. In contrast 
to the IVSGDR, which displays a complex underlying struc- 
ture with three overlapping components and its strength spreads 



Table 1 : Values of the neutron-skin thickness (ARp„) of Sn determined using 
various experimental methods, in comparison with the neutron-skin thickness 
deduced in the present work. 



Method 


Ref. 


Date 


AR„„ (fm) 


{p,p) 0.8 GeV 


[8] 


1979 


0.25 + 0.05 


(Q',a')IVGDR 120 MeV 




1994 


0.21+0.11 


antiproton absorption 




2001 


0.19 + 0.09 


Cne,t) IVSGDRh-AGDR 


[10] 


2004 


0.27 + 0.07 


pygmy dipole resonance 


U5,i3J 


2007 


0.19 + 0.05 


{p,p) 295 MeV 




2008 


0.185 + 0.05 


AGDR present result 




2013 


0.21 +0.05 



over a large energy interval, the AGDR represents a rather sim- 
ple charge-exchange mode (A^ = 1", AL - I). While previous 
analyses were based on sum rules, this work introduces an al- 
ternative self-consistent approach that could systematically be 
used not only for the AGDR, but also for other modes sensitive 
to the neutron skin. 

As a first test, we have used the self-consistent RHB plus 
proton-neutron RQRPA to calculate the Tq - 1 dipole excita- 
tions in a sequence of Sn isotopes. By using effective interac- 
tions with density-dependent meson-nucleon couplings in the 
particle-hole channel, and the pairing part of the Gogny inter- 
action DIS for the T=l pairing channel, it has been possible to 
reproduce the experimental results on the excitation energy of 
the AGDR relative to the isobaric analog state. We have also 
shown that the isotopic dependence of the energy difference be- 
tween the AGDR and IAS provides direct information on the 
evolution of neutron-skin thickness along the Sn isotopic chain. 
Very good results have been obtained in comparison with avail- 
able data on the neutron-skin thickness. The present analysis 
demonstrates that this quantity can be determined by measur- 
ing the excitation energies of the AGDR relative to IAS. 

The accuracy of the method has been tested in the example of 
'^"^Sn. By employing a set of effective interactions that span a 
broad range of values for the neutron-skin thickness (as a result 
of variation of the symmetry energy at saturation density), the 
size of the neutron skin has been determined from the AGDR 
energies relative to IAS. The result is in very good agreement 
with previously published experimental values. More extensive 
studies, in line with recent work on the electric dipole polariz- 
ability and neutron-skin thickness |7] that has employed fami- 
lies of non-relativistic and relativistic energy density function- 
als, would allow a further reduction of theoretical uncertainties. 
The successful test of the method based on the AGDR holds 
promise for determining the size of the neutron-skin of unsta- 
ble neutron-rich exotic nuclei, and this is reinforced by recent 
advances in the development of RIBs. 
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